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ABSTRACT 

We demonstrate that the chaotic orbits of Prometheus and Pandora are due to inter- 
actions associated with the 121:118 mean motion resonance. Differential precession 
splits this resonance into a quartet of components equally spaced in frequency. Libra- 
tion widths of the individual components exceed the splitting resulting in resonance 
overlap which causes the chaos. A single degree of freedom model captures the es- 
sential features of the chaotic dynamics. Mean motions of Prometheus and Pandora 
wander chaotically in zones of width 1.8 deg yr _1 and 3.1 deg yr _1 , respectively. 
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1 INTRODUCTION 

Goldreich and Rappaport (2003) (hereafter abbreviated as GR) showed that the motions 
of Prometheus and Pandora are chaotic. The chaos arises from their mutual gravitational 
interactions, which explains why their longitude discrepancies have comparable magnitudes 
and opposite signs (French et al. 2002). Numerical integrations that account for the full 
mutual interactions and Saturn's gravitational oblatenss yield a Lyapunov exponent of order 
0.3 yr -1 . Although the results reported by GR assumed satellite masses based on a nominal 
density of 0.63 g cm" 3 , the Lyapunov exponent is insensitive to the assumed density above a 
critical value of approximately 0.3 gem" 3 . 

GR's integrations also reproduce qualitative features of the discrepancies between the 
longitudes of the satellites derived from analysis of recent HST data and predictions based 
on orbits fit to Voyager images (French et al. 2002). Sudden changes in the mean motions of 
Prometheus and Pandora are a striking feature of the numerical integrations. These occur 
at intervals of 6.2 yr when the satellites' apses are anti-aligned. It is notable that the only 
clearly documented changes in the mean motions occurred around the time of the most 
recent apse anti-alignment (cf. GR). 

The plan of this paper is as follows. In §2 a quartet of 121:118 mean motion resonances is 

identified as the probable cause of the chaos. 1 Then we describe two new programs in which 

interactions between the satellites are limited to those due to this quartet. The simpler of 
1 Differential apsidal precession splits each mean motion resonance into a multiplet of closely spaced 

components. 
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these reduces the resonant dynamics to a system with one degree of freedom. Results from 
these programs are compared in §3 with those obtained from integrations that account for 
the full gravitational interactions. §4 is devoted to a discussion of the width of the chaotic 
zone. 



2 ORIGIN OF CHAOS 

2.1 Resonant Quartet 

A systematic search for j : j — k mean motion resonances with k < 4 turned up j = 121, k = 
3. 

Following Murray and Dermott (2001), we write the disturbing function for the action 
of Pandora on Prometheus as 

^ Gin' , . 

K=—R D , 1 
a' 

and that for the action of Prometheus on Pandora as 

W = ^Ro • (2) 

a 

Here m and a denote mass and semi-major axis, and G is the gravitational constant. 2 To 
lowest order in the eccentricities, the terms in the literal expansion of the disturbing function 
associated with a k = 3 resonance take the form 

R D = e 3 f 8 2 cos [121 A' - 118 A - 3zu] + 



2 Primed and unprimcd symbols refer to Prometheus and Pandora, respectively. 
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eV / 83 cos [121A' - 118A - 2zu - w'\ + 
ee' 2 / 84 cos [121A' - 118A -w- 2zu'] + 
ee' 3 / 85 cos [121A' - 118A - 3tu'] , (3) 

where e, A, and w stand for eccentricity, mean longitude, and argument of periapse. The 
fsn are expressed in terms of Laplace coefficients evaluated at a = a/a'. 3 



fs2 = ^ | (-26j + 30j 2 - 8j 3 ) + (-9 + 27 j - 12j 



2 x db[%(a) 

" 1 a — '- h 

da 



/ss = ^ I (" 9 + 31j " 30j2 + 8jS ) + (9 - 25 j + 12j 2 ) a 



16 T •' •' •' y 1 - v - •' •' ' r/n 

- 5 + ^ a2 ^ i + « 3 ^ i ^ 



1 32 J + 30j 2 -8j 3 )'^ 2 ^ ^ ' ' ^ 2)(a) 

(4 - 6j > 2 ^^ - a 3^2^ I (6) 



/ 8 4 = - ] (8 - 32j + 30j 2 - 8j 3 ) b^ /2 2 \a) + (-8 + 23 j - 12j 2 ) a ^ 

,rf 2 ^ 2) (a) _ 3 d 3 b[y 2 2) (a 
da 2 da 3 



db ( L 2) (a) 



,d 2 b^ 2 3 \a) d 3 b^ 2 3 \a)} 

Tables 1 and 2 list values for the parameters used in this paper. Satellite masses are given 

as fractions of Saturn's mass based on an assumed density of 0.63 gem -3 . Initial values for 
3 In the following, j should be viewed as a shorthand for 121. 
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mean longitudes, apsidal angles, mean motions, and eccentricities are based on orbits fit 
to Voyager images by Jacobson (2001). Precession rates are calculated from the Saturnian 
gravitational field (Campbell and Anderson 1989). 

Table 1 . Masses, Initial Mean Longitudes, & Mean Motions. 



Satellite 


m/M 


Mean Longitude (°) 


Mean Motion (°/s) 


Prometheus 
Pandora 


5.80 x 10~ 10 
3.43 x 10~ 10 


188.53815 
82.14727 


6.797331 x 10~ 3 
6.629506 x 10" 3 


Table 2. Eccentricities, Initial Apsidal Angles, & Precession Rates. 


Satellite 


Eccentricity 


Apsidal Angle (°) 


Precession Rate (°/s) 


Prometheus 
Pandora 


2.29 x 10~ 3 
4.37 x 10~ 3 


212.85385 
68.22910 


3.1911 x 10~ 5 
3.0082 x 10~ 5 



Rates of change of the arguments, corresponding periods, and coefficients for the four 
terms in equation Q are given in Table 3. 

Table 3 . Resonance Arguments, Rates of Change, Periods, and Coefficients. 
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Argument 


Rate of Change 


(°/s) 

V / / 


Period (year) 


Coefficient 


12LV- 118A 


— 3w 


-1.058 x 10" 


-5 


1.078 


e 3 / 82 = -0.001 


121A'- 118A 


— 2t*7 — ccj' 


-0.875 x 10- 


-5 


1.303 


e 2 e'fs3 = 0.006 


121A'- 118A 


— zu — 2m' 


-0.692 x 10- 


-5 


1.648 


ee ,2 / 8 4 = -0.01 


121A'- 118A 


-3m' 


-0.509 x 10- 


-5 


2.239 


e' 3 / 85 = 0.00 7 



2.2 Numerical Integrations 

To demonstrate that the quartet of 121 : 118 resonances is the cause of chaos in the 
Prometheus Pandora system, we develop two new programs to integrate the satellites' equa- 
tions of motion. Interactions between the satellites are restricted to the four resonant inter- 
action terms in the Fourier expansion of the disturbing function R D . Each program evolves 
propagates the satellites' orbital elements rather than their cartesian coordinates and veloc- 
ities as is done by the "old program" FSHEP used in GR. 

We adopt epicyclic elliptic elements instead of the more standard osculating elliptic el- 
ements since, unlike the latter, they do not require short period terms to describe elliptic 
orbits around oblate planets (cf. Borderies-Rappaport and Longaretti 1994; henceforth, re- 
ferred to as BRL.). BRL derive a modified version of Gauss' equations for the elements a e , 
e e , m e = u e + Q e , and A e = m e + M e . 4 From these, it is a straightforward exercise to derive 
a restricted version of Lagrange's equations that is valid in the planar case. We work with 

a simplified set appropriate for low eccentricity orbits about a modestly oblate planet. The 
4 Hcreafter we drop the subscript e. 
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equations read: 



~dt 



n 



da 
~dt 



2 dK 



(9) 



eft 



(10) 



where 



eft 



GM 


1 + 


3 




a 3 


2 


V a 


GM 


1 - 


3 


'Rp 


a 3 




2 


K a 



■h 



i on 



Ka 2 e dw 



15 (R V \ A 35 /i2p 
— \ J\ + 

a / 



a / 



16 V a 
175 /H 



16 



J fi 



:n) 



(12) 
(13) 



2.3 New Programs 

In the planar approximation, the Prometheus-Pandora system has four degrees of freedom 
and preserves two integrals, total energy and total angular momentum. Thus each phase 
space trajectory lies on a six dimensional hypersurface embedded in the eight dimensional 
phase space. 

FSHEPRES integrates the four, first-order equations (fH|)-([TTj) for each satellite. Thus 
it differs from FSHEP mainly because it limits the interactions between the satellites to 
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resonant terms. 5 Other minor differences arise because FSHEPRES integrates a simplified 
set of Lagrange's equations. In particular, the conservation laws are only approximately 
satisfied. 

FSHEPSIM integrates only the first two equations (jSJ) and Q for each satellite. This 
drastic simplification is reasonable because, as a consequence of the rapid differential pre- 
cession caused by Saturn's oblateness, interactions between the satellites produce negligible 
effects on their apsidal angles and orbital eccentricities (GR). A further simplification arises 
because conservation of energy implies 6 

m da m' da' , . , 

cr dt a dt 

Thus the resonant dynamics of the Prometheus-Pandora system reduces to a single degree 
of freedom system. It proves convenient to define the variable 

i> = 121A' - 118A , (15) 

so that Rd is expressed as 

4 

R D = J2 c kOos(^ + 5 k ) , (16) 

k=l 

with each 5 k = constant. The evolution of ip is governed by 

4 

£C*sm(^ + <f fc ) (17) 

k=l 
4 

= 3x (12in') 2 ^-[l + a(m'/m)]J2C k sm(^ + 5 k ) , (18) 

M k=i 

5 We view as an unimportant difference the use of orbital elements by FSHEPRES and cartesian positions 

and velocities by FSHEP. 

6 In the following equation we ignore the interaction energy which is only significant near conjunctions. 
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dt 2 



where in writing the second form of equation (JTHJ), we have applied the mean motion reso- 
nance relation Q'/Q « 118/121 and emphasized the contribution from the lighter body, ml . 
Individual mean longitudes follow from the relations 

() -aim 1 /m)jj(t) + 118[A(0) + A(0)t] + 121a(m'/m)[\'(0) + A'(0)t] ( , 
^ ~ 121[1 + a{m'/m)} ' ^ ' 

m + 121a(m'/m)[A'(0) + A'(0)t] + 118[A(0) + A(0)t] 
U 121[l + a(m'/m)] ' 1 J 

Although we are left with a system described by a single degree of freedom, the absence of 
an energy integral still allows for chaos. 



3 COMPARISON OF RESULTS 

In this section we compare result obtained using FSHEPRES and FSHEPSIM with those 
obtained with FSHEP. As in GR, all our simulations are initialized with orbital elements 
for Prometheus and Pandora taken from Jacobson's ephemerides at epoch 1981 August 23 
04:02:12 UTC. Comparisons among similar calculations done with each of the three programs 
are presented in Figures 1-6. As a consequence of chaos, qualitative similarities are the best 
that can be expected. These are apparent in each set of figures. However, there is a hint that 
FSHEPSIM yields a slightly smaller Lyapunov exponent than either FSHEP or FSHEPRES. 
The similarity between the 20 year runs of longitude variations displayed Figures 1 and 

2, while consistent with a Lyapunov exponent of approximately O.Syr" 1 as shown in Figure 

3, probably also reflects the fact that at the Voyager epoch the mean motions of Prometheus 
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and Pandora were close to their respective maximum and minimum. This accounts for the 
negative values of the rates of each resonant argument quoted in Table 3. 

Figure 4 shows that over 3000 years the net variation of 121 A' — 118A — 2m' — m is 
much smaller than that of the other phases. Together with the constraint imposed by the 
conservation of energy on relative variations of n and n', this implies that over this time 
interval the average mean value of n was smaller than its initial value by about 0.67 deg yr" 1 
and that of n! was larger by about 1.14 deg yr _1 . The relatively small value of 121n' — 
118n — 2m' — w has a plausible dynamical explanation. The term with this phase rate is the 
one with the largest amplitude. Moreover, the amplitudes of terms with neighboring phase 
rates are about half as large and have opposite signs to that of the dominant term, whereas 
the amplitude of the term with the slowest phase rate is much smaller. Support for this 
explanation is provided by observing that interchanging the values of e and e' results in the 
phase 121A' — 118A — w' — Izu assuming the special status of being the one with the smallest 
net variation. 

Figures 5 and 6 display longitude variations over 3000 years relative to the longitude 
based on the average mean motion over this interval. These are seen to be bounded by 
±180 degrees. 

4 Discussion 

A closer examination of the one-degree of freedom model developed for FSHEP provides 
additional insight regarding chaos in the Prometheus-Pandora system. 
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Overlapping resonances are known to produce chaos. Frequencies of individual mem- 
bers of the resonant quartet are spaced by w — m' ~ 1.0 radyr -1 . This is smaller than 
the half widths of the individual resonance components. 7 Half widths computed from 
equation (|T5|) and the data in Tables 1-3 are, in order of increasing resonance frequency, 
1.5, 3.7, 5.1, 4.1 radyr" 1 . 

Figure 7 shows surfaces of section based on data from 3000 year integrations using 
FSHEPSIM. A single point with coordinates if) — m — 2m', if) — m — 2m' is plotted each 
time the apses align (every 6.2 yr when m — m' = modulo 2tt). 8 Nominal values for the 
satellites' masses were used for the upper panel. The scattering of points over an area in the 
phase plane is a signature of chaos. The balance in the number of points above and below the 
horizontal axis and the overall vertical width of their distribution are a consequence of the 
dominance of the resonance component with phase 121A' — 118A — m — 2m'. Satellite masses 
were reduced by a factor 10 below their nominal values to obtain the integrable example 
whose surface of section is shown in the lower panel. 

Variations of n and ml are related to those of if> by 



dn 
~dt 



1 d^_ 

121 [1 + a{m'/m)]~dt 



(21) 



dn' 
~dt 



118[1 + a(m'/m)] dt 



—a(m'/m) dip 



(22) 



Thus the total width in if) corresponds to full width variations An m 1.8degyr 



1 and 



7 The half width is the maximum angular velocity achieved during motion on the separatrix. 
8 We chose apse alignment to minimize the effects of the interaction energy. 
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An' 3.1 deg yr 1 . 
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7 Figure Captions 



FIGURE 1 : Prometheus longitude in degrees from numerical integration as a function of 
time over 20 years. A drift based on the initial mean motion has been subtracted from the 
longitudes. Dashed lines indicate the times of periapsis antialignment. Results shown in the 
top, middle, and bottom panels were obtained with the programs FSHEP, FSHEPRES, and 
FSHEPSIM. 

FIGURE 2 : Pandora longitude in degrees from numerical integration as a function of time 
over 20 years. A drift based on the initial mean motion has been subtracted from the 
longitudes. Dashed lines indicate the times of periapsis antialignment. Results shown in the 
top, middle, and bottom panels were obtained with the programs FSHEP, FSHEPRES, and 
FSHEPSIM. 

FIGURE 3 : Lyapunov exponent in yr" 1 for the Prometheus-Pandora system over a period 
of 3 x 10 3 years. The results shown in the top, middle, and bottom panels were obtained 
with the programs FSHEP, FSHEPRES, and FSHEPSIM. 

FIGURE 4 : Phases ip + 5k, k — 1, . . . , 4, in radians along the solution. Results shown in the 
top, middle, and bottom panels were obtained with the programs FSHEP, FSHEPRES, and 
FSHEPSIM. 

FIGURE 5 : Prometheus longitude in degrees from numerical integration as a function of 
time over 3000 years. A drift based on the mean motion averaged over 3000 years has been 
subtracted from the longitude. Results shown in the top, middle, and bottom panels were 
obtained with the programs FSHEP, FSHEPRES, and FSHEPSIM. 

FIGURE 6 : Pandora longitude in degrees from numerical integration as a function of time 
over a 3000 years. A drift based on the mean motion averaged over 3000 years has been 
subtracted from the longitude. Results shown in the top, middle, and bottom panels were 
obtained with the programs FSHEP, FSHEPRES, and FSHEPSIM. 

FIGURE 7 : Surfaces of section obtained by plotting (ip— w— 2m', ip—vo—2vo') at each time of 
periapsis alignment over 3,000 years. Units are radians and radians per year. Computations 
were made with FSHEPSIM, for the top panel with the nominal value of 0.63 gem -3 for the 
satellite density. For the bottom panel, the density was reduced by a factor of 10 in order 
to obtain an integrable example to contrast with the chaotic one shown above. 
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